Ferromagnetism of Weakly-Interacting Electrons in Disordered Systems 
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It was realized two decades ago tiiat tiie two-dimensional diffusive Fermi liquid phase is unstable against arbi- 
trarily weak electron-electron interactions. Recently, using the nonlinear sigma model developed by Finkelstein, 
several authors have shown that the instablity leads to a ferromagnetic state. In this paper, we consider diffusing 
electrons interacting through a ferromagnetic exchange interaction. Using the Hartree-Fock approximation to 
directly calculate the electron self energy, we find that the total energy is minimized by a finite ferromagnetic 
moment for arbitrarily weak interactions in two dimensions and for interaction strengths exceeding a critical 
value proportional to the conductivity in three dimensions. We discuss the relation between our results and 
previous ones. 
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I. INTRODUCTION 

Since the discovery of Anderson localization, it has been 
well known that disorder can profoundly change the electronic 
properties of solids. For non-interacting electronic systems, 
the celebrated scaling theory of Abrahams et.ali shows that 
electrons are always localized in one and two dimensions, 
while in three dimensions, there exists a metal-insulator tran- 
sition at a critical value of the disorder strength. On the other 
hand, when interactions between electrons are taken into ac- 
count, the results are much less clear In a disordered sys- 
tem, at length scales beyond the mean free path, electrons 
move diffusively. Since diffusion is slow, electron-electron 
interactions are consequently enhanced. First, electrons' abil- 
ity to screen the long-ranged Coulomb interaction is ham- 
pered. Second, electrons stay in each other's region longer 
and hence interact with each other more strongly^ As shown 
by Finkelstein^ and Castellani et al'^, the two-dimensional dif- 
fusive Fermi liquid phase is, as a result, unstable. However, it 
is not clear what is the stable fixed point controlling the low- 
energy behavior of weakly-interacting electrons in the pres- 
ence of disorder. 

Using the nonlinear sigma model developed by 
Finkelstein^, several authors^-*-^ have shown that the en- 
hancement of electron-electron interactions due to disorder 
leads to the formation of ferromagnetic moment. The 
non-linear sigma model takes as its starting point a repulsive 
density-density interaction which is customarily split into 
small-angle and large angle scattering with coupling constants 
Fi and r2^. The latter coupling drives ferromagnetism. In 
this paper, we consider itinerant electrons interacting through 
a ferromagnetic exchange iJjnt = — / — r')S(r) • S(r'). 
As in the case of a density-densty interaction, this is an 
interaction which, in the clean, disorder-free case, does not 
lead to ferromagnetism for small J; if it occurs at all, it can 
only be at large J. This is because a kinetic energy penalty 
must be paid when the up- and down-spin Fermi surfaces 
are split in order generate a net ferromagnetic moment. 
The kinetic energy penalty for splitting the Fermi surfaces 
by (5/i is ^ vSji^, where v is the density of states at the 
Fermi surface. On the other hand, the energy gain for non- 
interacting electrons is ^ Ji/^Sfi^. Hence, ferromagnetism 
can only occur for J > -i, the Stoner criterion. However, 



as we shall see in the next section, in a two-dimensional 
disordered system, the interaction energy gain varies instead 
as ~ Jv^ ln(r5/x); consequently, it is always favorable 
to develop a ferromagnetic moment. (This is similar to 
what occurs in a quantum dot**.) We show this by directly 
calculating the disorder-averaged electron self-energy, using 
the trick of ref. 9, in which it is related to the diffusive pole 
in the density-density correlation function. This method is 
physically transparent since one can directly see how the 
seemingly innocent diffusive pole leads to a logarithmic 
correction to the self energy. We also avoid the need for the 
replica trick used in the sigma model treatment. For a three 
dimensional system, our calculation shows that ferromagnetic 
order develops when J exceeds a critical value which is 
proportional to the conductivity of the system. We compare 
these results with previous ones. 



II. THE MODEL 

We consider the following Hamiltonian for a weakly- 
interacting electronic system in the presence of a random dis- 
tribution of impurities: 



H — Ho + -ffir 



where 



Hq — ^ E„ 



(1) 



(2) 



is the noninteracting Hamiltonian, which includes the kinetic 
energy of electrons and the random potential. £',„, (pm{r) 
are respectively the m'^ eigenvalue and eigenfunction of Hq, 
which we assume to be spin-independent. We will call Em the 
'kinetic' energy of the electron to distinguish it from the inter- 
action energy associated with Hint, even though Em includes 
the effects of the random potential. Cma is the annihilation op- 
erator for an electron of spin a —1, | in the eigenstate Lp„i{r)- 
The annihilation operator tpa (r) for an electron of spin a at r 
is: 



(3) 
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As mentioned in the introduction, the interaction Hamihonian 
i/iiit considered here is: 

i?int ^^Jd'^r d''r' J(r - r')S(r) • S(r') 

+ 2{S+{r)S-{r') + S-{r)S+{r'))) (4) 



down-spin electrons at the Fermi level: M = v Sji, where v 
is the single-electron density of states. Now, we find that the 
disorder-averaged variational interaction energy is: 



(vPol i/int |*o> = - y d'^vd^v' J(r ~ r') {S^-{r)) (5-(r')) 
j dEdE' j d'^rd'^r' J{r-r')F{E,E',r,r') NeN, 



Here J(r — r') > is the exchange interaction which is as- 
sumed to be weak, ferromagnetic, and short-ranged. S(r) = where 
ipl^{r)(Ja0ip(3{r) is the spin density operator, and a are the 
Pauli matrices. 



E' 

(10) 



(S^ir)) - J2^*„,{r)ipmir) {N^^ - N^^) (11) 



III. HARTREE AND FOCK ENERGIES 

As a warmup, let us first consider a clean system, for which 
the ifm{r) are plane waves. We look for ferromagnetism in 
this model by assuming taking a trial wavefunction l^'o) clean 
with a ferromagnetic moment: 



(5) 



and minimizing ('J'ol H j'J'o) with respect to the magnetiza- 
tion M ^ V Sji, where v is the single-electron density of states 
and 2 (5/i = {k'^-\ ~ ^fx) /2™ is the difference in the kinetic 
energies of up- and down-spin electrons at the Fermi energy. 

In evaluating the expectation value of the interaction en- 
ergy, we find two terms: 

('J'ol ^fintl^o) --^(0)1/2 + 

^ J(k - k') {^-NuNk' - \SISI^ (6) 



k,k' 



where iV/c, S*^ = 9{kFi—k)±6{kFi — k). They are essentially 
the Hartree and Fock contributions to the energy. 
If we take for simplicity J(q) = J, then 



Z o 



(7) 



Ferromagnetism will only occur if this interaction energy is 
larger than the kinetic energy cost 



(-ffkin) M ~ (-f^kin)o = ^ 



1 



(8) 



i.e. only for 3 J/2 > v. (The magnetization is stabilized at a 
finite value by the energy dependence of the density-of states, 
which we have neglected here.) 

Consider, instead, the dirty case. We take 



I* 



OMirty^ n n cI^ci,JO) (9) 



where En„ = Ep. Again, the magnetization M is propor- 
tional to the difference in the 'kinetic' energies of up- and 



{Nmtl) = (cLta^^Ta) = ± M - m) (12) 

Ne is the number of electrons in the state at energy E, namely 

TVs = 2 for £; < ~Sfi, Ne ^ 1 for ~Sfi < E < Sfi, and 
Ne = for E > S^i. The function F{E, E'; r, r') is defined 
as 

FiE, E'- r, r') ^Y,5{E - E^)S{E' - E^) x 



^*miT)^m{r')^pl{r')^n{r) (13) 

The upper bar indicates average over disorder. The first term 
in is similar in behavior to the clean case, and will be 
neglected in the following. The second term is more interest- 
ing because the function F is closely related to the density- 
density correlation function. The crucial difference with the 
clean case is that the latter has a diffusive form, as we now 
describe. 

We assume that the density propagates diffusively, so that 



n(q,w) = i / dtd''x([p(x,t),p(0,0)])e 

Dq^ 

= V ij 



(14) 



where, again, v is the density of states at the Fermi surface, 
and D is the diffusion constant. Then 



n*(q,c^)=-i / did'*x([p(0,0),p(x,t)])e' 
Jo 

ro 

= -i dtd''x([p(x,<),p(0,0)])e~'i-^+''^* (15) 

^ — oo 

On the second line, we have used the fact that after averag- 
ing over disorder, the density-density correlation function is 
translationally invariant. This assumption is valid in the high 
conductivity regime of interest here. From this, we can ob- 
tain the Fourier transform of the density-density correlation 
function. 



dtd'^xe 



([p(x,t),p(0,0)]) (16) 
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by observing that 



1 



Aiq,u;) = -(n(q,c^)-n*(q,c^)) (17) 
= ilm(n(q,^)) 

We can relate this to our function F by noting that 

{[p{rt),p{r',Q)]) ^Y.^*rn{r)v,n{r')^Ur>n{r) 

x{N^-N„) (18) 

So we get 

poo pfJ. 

A{q,uj) ^ dE dE' F{E,E';r,r') 

J ^ J — OO 

X e"i^''"''H{E - E' - uj) 
= ^F(si,uj) (19) 

Where, on the last line, i^(q, is the Fourier transform of 
F(E^ E'; r, r'), and lu is the energy difference E — E'. /i is 
the chemical potential. So we get 



and the second term in (I10> becomes 



(20) 



(*o|i?i„t|*o> 



/ dE r dE' ^ '^"'^ 



Dq^ 



(27r)d {E - E')^ + {Dq^Y 



J(q) (21) 



The . . . refers to the first term in (llOt . which is smaller because 
it does not have any infrared singularities in its integral. In the 
preceding derivation, we performed the disorder average of F, 
by comparing it to the density-density correlation function, 
thereby avoiding the replica trick normally used in the field 
theoretical treatment. 



IV. FERROMAGNETIC GROUND STATE 

We now minimize {H) with respect to the magnetization 
M. We define the effective interaction V{E — E') between 
electrons in states at energies E, E' by 

(*o|-ffint|*o) = / dE dE' V{E - E') NeNe' (22) 

J E>E' 

Evaluating the momentum integral in \2l\ . we obtain 



V{E-E') = 



vJ 



In 



{E-E'f + A^ 



{E- E- 



•l\2 



(23) 



Here, we have ignored the momentum dependence of J(q), 
and as usual, A is the ultraviolet cutoff which can be taken to 
be 1/t, where r is the elastic scattering time. 



Thus, the interaction energy is: 



dE / dE' V{E ~ E')-2-2 

-A J-A 

/Sfj, p-Sn 
dE / dE' V{E -E')-2-l 

/Sp pE 
dE 1 dE' V{E -E')-l-l (24) 

Subtracting the interaction energy in the paramagnetic state, 

M = 0, 

^0 pE 

(i/int>o = / dE dE' V{E -E')-2-2 (25) 

J-A J-A 



we have 



iiit/o 



dE I dE'V{E~E') 



vJ 



5}.L 



dE 



dE' In 



{E - E'f + A2 



vJ 
8^ 



dE{E + 5n) ^1 + In 



{E-E'f 
{E + Sfif + A2 



v.] 

'S,ttD 



{E + 5pf 



A 



64i^ + 44i^ln — 

2(3/i 



(26) 



Including the 'kinetic energy', we have 



vJ 



A \ 



(27) 



Even at small J, ferromagnetic order occurs because the 
total energy is minimal at a finite value of Sji and, therefore, 
of M: 



2 



(28) 



So we have shown that the diffusive fermi liquid is always 
unstable against the ferromagnetic order in two dimensions. 

In three dimensions, one can do the analogous calculations, 
leading to 



{H) - Eo 



2 JM^^,2 



4Jj/ 



5/2 



(29) 

In the first term, we have included the Hartree contribution 
because the Fock part is non-divergent and can be viewed as 
a correction. Hence, the critical coupling Jc above which fer- 
romagnetism occurs is modified from the Stoner value: 



(30) 



where £ = vpris the elastic mean-free path and we have used 
D = VpT/ 3 in three dimensions. 
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V. DISCUSSION 

In this paper, we use a straightforward Fock approximation 
to show that in two dimensions an arbitrarily weak ferromag- 
netic interaction causes an instabiHty of the diffusive Fermi 
hquid. In three dimensions, this instability occurs only when 
the interaction strength exceeds a critical value proportional 
to the conductivity. We worked directly with the electronic 
Hamiltonian, without recourse to Finkelstein's non-linear a 
model, which was used at an intermediate stage of previous 
calculations^-^-^. Our results are consistent with these earlier 
ones. The model and method of this paper has the advantage 
of physical transparency: it shows the direct connection be- 
tween diffusive motion and the ferromagnetic instability. 

By using the method in ref. 9, we avoided using the replica 
trick, and saw the crucial feature of diffusive motion: it im- 
plies a spatial correlation between single-particle states which 
are nearby in energy. As a result of this correlation, the effec- 
tive interaction strength diverges as the energy of two elec- 



trons approach each other This strongly favors ferromag- 
netism since electrons of the same spin cannot occupy the 
same energy level. This strong dependence on the energy sep- 
aration of the electrons is completely absent in a clean system. 

We note that once ferromagnetic order develops, the system 
goes into the universality class of electrons in a Zeeman fieldi, 
which flows towards an insulating state in two dimensions. 
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